Discontinuous transitions in double exchange materials 
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It is shown that the double exchange Hamiltonian, with weak antiferromagnetic interactions, has a 
rich variety of first order transitions between phases with different electronic densities and/or mag- 
netizations. For band fillings in the range 0.3 < x < 0.5, and at finite temperatures, a discontinuous 
PM-FM transition between phases with similar electronic densities but different magnetizations 
takes place. This sharp transition, which is not suppressed by electrostatic effects, and survives in 
the presence of an applied field, is consistent with the phenomenology of the doped manganites near 
the transition temperature. 
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I. INTRODUCTION. 

Doped manganites show many unusual features, the 
most striking being the colossal magnetoresistance 
(OMR) in the ferromagnetic (FM) phase |]-||. In addi- 
tion, the manganites have a rich phase diagram as func- 
tion of band filling, temperature and chemical composi- 
tion. The broad features of these phase diagrams can 
be understood in terms of the double exchange model 
(DEM) j§[|, although Jahn-Teller deformations [g and 
orbital degeneracy may also play a role 0]. A remark- 
able property of these compounds is the existence of in- 
homogeneities in the spin and charge distributions in a 
large range of dopings, compositions and temperatures 
|^| ^0[ . At band fillings where CMR effects are present, 
x ~ 0.2 — 0.5, these compounds can be broadly classified 
into those with a high Curie temperature and a metallic 
paramagnetic phase, and those with lower Curie temper- 
atures and an insulating magnetic phase [Tl|-[r^|. 

The DEM is a simplification of the FM Kondo lattice, 
where the FM coupling between core spins and conduc- 
tion electrons is due to Hund's rule. When this cou- 
pling is larger than the width of the conduction band, 
the model can be reduced to the double exchange model 
with weak inter-atomic antiferromagnetic (AFM) inter- 
actions. Early investigations juj showed a rich phase 
diagram, with AFM, canted and FM phases, depending 
on doping and the strength of the AFM couplings. More 
recent studies have shown that the competition between 
the double exchange and the AFM couplings leads to 
phase separation into AFM and FM regions, suppressing 
the existence of canted phases [jl5|— |l7]Jl9|l . In addition, 
the double exchange mechanism alone induces a change 
in the order of the FM transition, which becomes of first 
order, and leads to phase separation, at low dopings [^0| . 
Note, however, that a detailed study of the nature of the 
transition at finite temperatures is still lacking, despite 
its obvious relevance to the experiments. 

The purpose of this work is to investigate systemat- 
ically the phase diagram of the DEM with weak AFM 



interactions. We find, in addition to the previously dis- 
cussed transitions, a PM-FM first order transition near 
half filling, if the double exchange mechanism is suffi- 
ciently reduced by the AFM interactions. This transi- 
tion does not involve a significant change in electronic 
density, so that domain formation is not suppressed by 
electrostatic effects. 

The model is described in the next section, and the 
method of calculation is introduced in the following sec- 
tion. The main results are presented in section IV, and 
the main conclusions are discussed in section V. 



II. THE MODEL. 

We study a cubic lattice with one orbital per site. At 
each site there is also a classical spin. The coupling be- 
tween the conduction electron and this spin is assumed 
to be infinite, so that the electronic state with spin an- 
tiparallel to the core spin can be neglected. Finally, we 
include an AFM coupling between nearest neighbor core 
spins p8[. The Hamiltonian is: 



(hi) 



JafS Si ■ Sj 



(1) 



where S — 3/2 is the value of the spin of the core, 
Mn 3+ , and S stands for a unit vector oriented paral- 
lel to the core spin, which we assume to be classical. 
In the following, we will use Jaf = JafS 2 . Calcula- 
tions show that the quantum nature of the core spins 
does not induce significant effects |T7|. The function 



T(S t , Sj) = t [cos f cos ^ + sin ^ sin ^-e 1 ^* 
for the overlap of two spin 1/2 spinors oriented along the 
directions defined by Si and Sj, whose polar and az- 
imuthal angles are denoted by and <p, respectively. We 
study materials of composition Lai-^Af^MnOa, where M 
is a divalent ion, and x < 0.5. In this composition range, 
the probability of finding two carriers in neighboring sites 
(two contiguous Mn 4+ ions) is small, so that a carrier 



stands 



1 



in a given ion has all the e g orbitals in the next ions 
available. Then, the anisotropies associated to the dif- 
ferences between the two inequivalent e g orbitals should 
not play a major role. On the other hand, if x > 0.5, 
we expect a significant dependence of the hopping ele- 
ments on the occupancy of orbitals in the nearest ions. 
In this regime, the equivalence of the two e g orbitals in 
a cubic lattice can be broken, leading to orbital ordering 
|7|,pT| (see, however p^). We will show that the main 
features of the PM-FM phase transition, for x < 0.5, can 
be understood without including orbital ordering effects. 
Moreover, in this doping range, anisotropic manganitcs 
show similar features p3j-|26|] , which suggest the existence 
of a common description for the transition. We will also 
neglect the coupling to the lattice. As mentioned below, 
magnetic couplings suffice to describe a number of discon- 
tinuous transitions in the regime where CMR effects are 
observed. These transitions modify substantially the cou- 
pling between the conduction electrons and the magnetic 
excitations. Thus, they offer a simple explanation for 
the anomalous transport properties of these compounds. 
Couplings to additional modes, like optical or acousti- 
cal phonons [ p7[ , and dynamical Jahn- Teller distortions 
p8| will enhance further the tendency towards first order 
phase transitions discussed here. 

We consider that a detailed understanding of the role of 
the magnetic interactions is required before adding more 
complexity to the model. 



III. METHOD. 

At finite temperatures, the thermal disorder in the ori- 
entation of the core spins induces off-diagonal disorder 
in the dynamics of the conduction electrons. The calcu- 
lation of the partition function requires an average over 
core spin textures, weighted by a Boltzmann factor which 
depends on the energy of the conduction electrons propa- 
gating within each texture. We have simplified this calcu- 
lation by replacing the distribution of spin textures by the 
one induced by an effective field acting on the core spins, 
which is optimized so as to minimize the free energy. 
The electronic energy includes accurately the effects of 
the core spin disorder on the electrons. Our calculation 
is a mean field approximation to the thermal fluctuations 
of the core spins, retaining, however, the complexity of 
a system of electrons with off-diagonal disorder. This 
approximation can be justified by noting that the con- 
duction electrons induce long range interactions between 
the core spins, that always favor a FM ground state. In 
general, our method is well suited for problems of elec- 
trons interacting with classical fields. 

In more mathematical terms, we have used the varia- 
tional formulation of the Weiss Mean-Field method |2^] 
to compute the free energy of the system. We first trace- 
out the fermion operators in ([!]), thus obtaining the ef- 
fective Hamiltonian for the spins, 



n cS ({s}) = j af J2 s * ■ s i 

{hi) 



(2) 



k B TV / dEg{E-{S}) log 



— 

1 + e T B r 



where g(E; {S}), is the fermionic density of states and 
V the volume of the system. The Mean-Field procedure 
consists on comparing the system under study with a 
set of simpler reference models, whose Hamiltonian H.q 
depends on external parameters. We choose 

H = - h i ■ S i ■ ( 3 ) 

i 

The variational method follows from the inequality 



T < J-q + (H c ff — Wo)o • 



(4) 



where J-q is the free energy of the system with Hamilto- 
nian (||), and the expectation values (-)o are calculated 
with the Hamiltonian 7io- The mean-fields {h} are cho- 
sen to minimize the right-hand side of (f|). The calcula- 
tion of the right-hand side of (4), requires the average of 
the density of states (see Eq.(3)) on spin configurations 
straightforwardly generated according to the Boltzmann 
weight associated to the Hamiltonian Tig and tempera- 
ture T. The key point is that g{E\ {S}) can be numer- 
ically calculated on very large lattices without further 
approximations using the method of moments 1 3^] (com- 
plemented with an standard truncation procedure plf). 
We have extracted the spin-averaged density of states 
on a 64 x 64 x 64 lattice (for these sizes, we estimate 
that finite size effects are negligible). For simplicity on 
the analysis, we have restricted ourselves to four families 
of fields: uniform, hi = h, giving rise to FM ordered 
textures; hi = (—l) Zi h, originating A-AFM order, i.e., 
textures that are FM within planes and AFM between 
planes; hi = (—l) Xi+Vi h, producing C-AFM order, that 
is, textures that are FM within lines and AFM between 
lines; and staggered, hi = (—l) Xi+Vi+Zi h, which origi- 
nate G-AFM order, i.e., completely AFM textures. We 
have chosen fields of these kind since they produce the ex- 
pected kinds of order, although this is not a limitation of 
the method. Once the spin-averaged density of states is 
obtained, it is straightforward to obtain the values of the 
mean- field that minimize the right-hand size of Eq.(|]), 
and the corresponding value of the density of fermions. 
Expressing the right-hand side of Eq.(||) as a function 
of the magnetization (or staggered magnetizations), we 
obtain the Landau's expansion of the free energy on the 
order parameter. 

It is finally worth mentioning when our calculation and 
the Dynamical Mean Field Approximation |l^,|32| are ex- 
pected to yield the same results. It is clear that the 
key point is the calculation of the density of states in 
Eq.(f|). For this problem of classical variables, the dy- 
namical Mean-Field is known to yield the same density 
of states than the CPA approximation f33[. Under the 
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hypothesis of spatially uncorrelated fluctuations of the 
spins, which holds in any Mean-Field approximation, the 
CPA becomes exact on the Bethe lattice with large co- 
ordination number. However, one cannot conclude that 
with our calculation we would get the same results on the 
Bethe lattice, since one has still to specify the probability 
distribution for the spins to be used in the CP A calc ula- 

™ the 



tion of the average density of states. In Refs. |32 
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calculation is done by identifying an effective Heisenberg- 
like mean field, which becomes exact when the magneti- 
zation is very small. Then, the distribution of spin orien- 
tations is equivalent to the one generated by an effective 
magnetic field. In this limit, our ansatz should reproduce 
the calculations reported in |}2 19 1, when implemented in 
a Bethe lattice. 

In order to study first order transitions, one must con- 
sider solutions at finite magnetizations. Then, the opti- 
mal Boltzmann weights need not coincide with the effec- 
tive field ansatz made here. Detailed DMFA calculations 
for the double exchange model |34|,|35|], however, show 
that the differences between the optimal DMFA distri- 
bution and that obtained with an effective field are small 
throughout the entire range of magnetizations. Thus, 
the scheme used in this work includes the same physical 
processes as the DMFA, but it is also able to describe 
effects related to the topology of the three-dimensional 
lattice, like those associated to the Berry phase, which 
arises from the existence of closed loops. Furthermore, 
the present scheme allows us to study the relative stabil- 
ity of phases, like the A and C antiferromagnetic phases 
described below, which can only be defined in a cubic 
lattice. 



IV. RESULTS. 

The model, Eq. ([!]), contains two dimensionless param- 
eters, the doping x, and the ratio JafA- The range 
of values of x is < x < 1, and the Hamiltonian has 
electron- hole symmetry around x = 0.5. The zero tem- 
perature phase diagram, shown in Fig. [j], is calculated 
minimizing the effective Hamiltonian at fixed chemical 
potential and zero temperature (we take the limit of 
zero temperature in Eq.(j3|) obtaining the grand-canonical 
Hamiltonian), within the four Mean-Field ansatzs pre- 
viously defined. At zero JafA, only the ferromagnetic 
phase is found, and the system is stable at all composi- 
tions. When Jaf A 1S finite, there is a value of the chem- 
ical potential for which the empty system with a perfect 
G-type AFM spin ordering has the same grand- canonical 
energy that a system with a perfect FM spin ordering and 
a finite value of a;. At this value of the chemical poten- 
tial the system is unstable against phase separation [^o) , 
as shown in Fig. [l]. Notice that the phase-separation re- 
gion can never reach x — 0.5, due to the hole-particle 
symmetry. For larger values of Jaf a small region of 
A-type AFM is found for x ~ 0.25, and a much larger 



region of C-type AFM for x close to half-filling. Finally, 
a G-type AFM-region is eventually reached by further 
increasing JafA- However, this is not a saturated anti- 
ferromagnetic phase since the mean-field that minimizes 
the grand-canonical energy has a finite h/T when T tends 
to zero |3^| (notice that one cannot have a continuously 
varying value of x in a perfect AFM configuration) . 




FIG. 1. Calculated phase diagram at T = 0. The A-AFM 
phase has ferromagnetic alignment within planes, and antifer- 
romagnetic alignment between parallel planes. The C-AFM 
phase has ferromagnetic alignment along chains, and antifer- 
romagnetic alignment between neighboring chains. 

Let us now discuss the phase diagrams at non-zero 
temperatures for the different values of JafA shown in 
Fig. |^. For Jaf = 0, we obtain a maximum transition 
temperature of T — 400K for a width of the conduction 
band W = \2t 2eV, which is consistent with a den- 
sity of states of p(Ep) = 0.85 eV _1 calculated in J3$| 
for La 1 / 3 Ca 2 / 3 Mn03 (see also j39|). Note that the band- 
width calculated in this way is probably an overestimate, 
as it does not include renormalization effects due to lat- 
tice vibrations f40|] . There is some controversy regarding 
the value of Jaf- The reported value of Jaf for the un- 
doped compound, LaMn03, is Jaf ~ 0.58meV, so that 
Jaf ~ 0.005t J4l[], although calculations give higher val- 
ues [Q. In the doped compounds, there is an additional 
contribution of order Jaf ~ t 2 /C/ ox , where U cyi « 1 — 2eV 
is the level splitting induced by the intra-atomic Hund's 
coupling [§. Thus, J A f ~ O.Olt- 0.08t, although higher 
values have been suggested . 

Our results show four types of first order transitions: 

i) In pure DEM (Jaf = 0) the magnetic transition 
becomes discontinuous at sufficiently low densities, in 
agreement with the analysis presented in pcfl . The phase 
coexistence region shrinks to zero and the critical tem- 
peratures vanish as x goes to zero, as expected. 

ii) The competition between antiferromagnetism and 
ferromagnetism when Jaf leads to a discontinuous 
transition which prevents the formation of canted phases, 
as reported in plH~7|. This transition also takes place 
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at low dopings. 




FIG. 2. Transition temperatures as function of electronic 
density and strength of the AFM couplings. The dashed lines 
correspond to continuous transitions. Solid thick lines are 
drawn for first order transitions, and the stripes correspond to 
phase coexistence regions. The onset of first order transitions 
at x ~ 0.5 is JafA w 0.06. 

iii) At moderate to high dopings, the FM-PM tran- 
sition becomes discontinuous, if the AFM couplings are 
sufficiently large. The onset for first order transitions at 
x = 1/2 is JafA ~ 0.06. Unlike the previous two cases, 
this transition takes place between phases of similar elec- 
tronic density. First and second order transition lines are 
separated by tricritical points. 

iv) In an interval of JafA, which depends on the dop- 
ing level, we also find phase transitions between the PM 
and A- AFM and C-AFM phases, that are of second or- 
der (see Fig. ||). At low temperatures there appear FM, 
C-AFM, A- AFM, and G-AFM phases separated by first 
order transitions with its associated phase separation re- 
gions, as shown in Fig |l|. 

As we see, the DEM complemented with AFM superex- 
change interactions between the localized spins give rise 
to a very rich magnetic phase diagram that contains first 
and second order transitions between phases with differ- 
ent magnetic order. 

In order to set a common frame for comparison with 
with standard approximations fl4| , po|| , we note the free 
energy of the system is made up of an entropy term, due 
to the thermal fluctuations of the core spins, an almost 
temperature independent contribution from the electrons 
and another temperature independent term due to the 
direct AFM coupling between the core spins. For in- 
stance, in the PM-FM case, we can write: T = 3 JafM 2 + 
E e \ ec (M)—TS(M) where S(M) is the entropy of a spin in 
an effective magnetic field producing magnetization M. 
We can expand: S(M) = -(§M 2 + ^M 4 + ^M 6 + ...) 
and E clcc (M) = dM 2 + c 2 M 4 + c 3 M 6 + ... where c±, c 2 
and C3 are functions of the band filling, and c\ is always 
negative [c\ , C2 and C3 are obtained fitting the numerical 



results for E c \ cc ). If there is a continuous transition, the 
critical temperature is given by Tq = (2|ci| — 6Jaf)/3. 
The transition becomes discontinuous when the quartic 
term in T{M) is negative. This happens if c 2 < and 
T < 20/9|c 2 |. Thus, if J AF > N/3 - 10|c 2 |/9, and 
c 2 < 0, the transition becomes of first order. A tricritical 
point appear in the transient between first and second 
order transitions. 

The fact that c 2 < is due to the energetics of the 
electrons in the disordered spin background. In a fully 
polarized system, M = 1, the electrons propagate in a 
perfect lattice. If M = the spins are completely disor- 
dered, and our results reduce to those reported in HJ45[ . 



Standard approximations 1 14 20J to the phase diagram 
of the DEM use the virtual crystal approximation, in 
which the cubic density of states is scaled by the average 
value {T(Si, Sj)}, defined in Eq.(|l|). This approxima- 
tion suffices to describe the main features of the phase 
diagram when Jaf = 0, but overestimates the kinetic en- 
ergy of the electrons moving in the disordered spin back- 
ground. The effect is more pronounced near half filling, 
when the electronic contribution is the largest, and c 2 is 
positive on the virtual crystal scheme. As our calculation 
takes fully into account the propagation of the electrons 
in a disordered environment, we think that the existence 
of a first order PM-FM transition when Tc is suppressed 
is a robust feature of the model. 
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FIG. 3. Transition temperature as function of the value of 
JafA f° r concentrations x = 1/3, x — 3/8, and x = 1/2. 
The dashed (solid) lines correspond to continuous (first or- 
der) transitions and a circle has been plotted at the tricrit- 
ical PM-FM point. In the top panel we sketch experimen- 
tal results from Ref. ^] where x = 1/3. The compound 
(Lai_ H Pr a ) 5 / 8 Ca 3 / g Mn03 studied in Ref. Q has x — 3/8. 
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At zero temperature, our calculation leads to a richer 
phase diagram to that calculated within the Dynamical 
Mean Field Approximation (4(|. As mentioned in the 
preceding section, our method coincides with this ap- 
proximation when implemented in a Bethe lattice. The 
topology of a cubic lattice allows for the possibility of 
A-AFM and C-AFM phases. 

We have developed an exact Monte Carlo algorithm to 
study the DEM. This approach is based in a Path Inte- 
gral formulation that allows to simulate on lattices much 
larger than in an usual Hamiltonian formulation. Full de- 
tails of the method will be given elsewhere . The first 
data of the Monte Carlo computation confirm the robust- 
ness of the present results. Simulations in the parame- 
ter region depicted in Fig. || show a very clear evidence 
for a first order transition in lattices up to 12 x 12 x 12 
sizes. In Fig. ^ we show data on a L = 8 lattice at 
half-filling at several temperatures. Note the large re- 
gion of metaestability marked by the vertical lines. It is 
also clear that fluctuations lower the transition temper- 
atures from their mean-field values, as it also happens 
in the three-dimensional Heisenberg model |Q. In ad- 
dition, we find a helicoidal spin structure at sufficiently 
low temperatures, which replaces, partially, the A-AFM 
and C-AFM phases discussed earlier. 
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FIG. 4. Monte Carlo results for the squared magnetization 
(bottom), and the k = (27r/8, 0, 0) squared Fourier compo- 
nent of the magnetization (top) in 8 x 8 x 8 lattices, as function 
of Jaf A, for x = 1/2 at different temperatures. 

Turning again to the Mean-Field approach, let us re- 
call that while a continuous transition is changed into a 
smooth crossover in an applied field, a first order tran- 
sition survives until a critical field is reached. The tran- 
sition takes place between two phases with finite, but 
different, magnetization, in a similar way to the liquid- 
vapor transition. The PM-FM line of first order transi- 
tions for dopings close to x = 0.5 ends in a critical point, 
(T C ,H C ). For Jaf = 0.08i, the critical field varies from 
H c = 0.00075t w 2.2T at x = 0.5 to H c = 0.0002< w 0.6T 
at x = 0.3, while T c w Tq and Tq is the Curie tempera- 
ture at zero field, shown in Fig. ||. 



V. CONCLUSIONS. 

We have shown that the phase diagram of double ex- 
change systems is richer than previously anticipated, and 
differs substantially from that of more conventional itin- 
erant ferromagnets. We have described first order tran- 
sitions which are either intrinsic to the double exchange 
mechanism, or driven by the competition between it and 
AFM couplings. In particular, we find that, in the doping 
range relevant for CMR effects, AFM interactions of rea- 
sonable magnitude change the PM-FM transition from 
continuous to first order. The existence of such a tran- 
sition has been argued, on phenomenological grounds, in 
order to explain the observed data in a variety (but not 
all) of doped manganitcs in the filling range x ~ 0.3 — 0.5 
[]49|,p0|| . The generic phase diagram that we obtain is 
consistent with a number of observations: 

i) Materials with a high transition temperature (low 
AFM couplings) have a continuous PM-FM transition, 
with no evidence for inhomogeneities or hysteretic effects. 
The paramagnetic state shows metallic behavior. 

ii) The PM-FM transition in materials with low transi- 
tion temperature (significant AFM couplings) is discon- 
tinuous. Near Tq inhomogeneities and hysteretic behav- 
ior are observed. The transport properties in the param- 
agnetic phase are anomalous. 

iii) Substitution of a trivalent rare earth for an- 
other one with smaller ionic radius (i.e., compositional 
changes that do not modify the doping level) diminishes 
the Mn — O — Mn bond angle, reducing the conduction 
bandwidth, W = 12t Assuming that the AFM 
coupling, Jaf, does not change significantly, the ratio 
Jaf A increases; therefore, the doping level y in series of 
the type (REi_ a RE y )i_ x AE x Mn03 might be traded by 
Jaf A- The top panel of Fig. |^ shows the experimental 
magnetic phase diagram of (Lai_ a Tb y ) 2 /3 Ca^ M11O3, 
as taken from Ref . [^3| . We note the similarities with the 
phase diagrams of the DEM in the plane (Jaf A, T/t) at 
fixed x. The phases A-AFM and C-AFM at intermediate 
Jaf A could become spin glass like phases in presence of 
disorder. 

iv) The first order PM-FM transition reported here 
survives in the presence of an applied field. A critical 
field is required to suppress it (hysteretic effects in an 
applied field have been reported in (52)). 
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